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ABSTRACT 

We build a matrix model of a chiral [SU{N)]^ gauge theory (SQCD5 decon- 
structed down to 4D) using random unitary matrices to model chiral bifunda- 
mental fields (N, N) (without (N,N)). We verify the duality by matching the 
loop equation of the matrix model to the anomaly equations of the gauge theory. 
Then we evaluate the matrix model's free energy and use it to derive the effective 
superpotential for the gaugino condensates. 
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1. Introduction 



Three years ago, Robbert Dijkgraaf and Cumrun Vafa discovered a peculiar duality between 
the A/" = 1 supersymmetric gauge theories in 4D and the bosonic matrix models without 
any spacetime at all. At first [1, 2], they showed that the effective superpotential for the 
gaugino condensates and the abelian gauge couplings follow from planar diagrams in the 
gauge theory (see also [3, 4]), and then they argued [5] that the very same perturbative 
series also gives the free energy of a bosonic random matrix model whose action is similar 
to the tree-level superpotential of the gauge theory. Shortly afterwards, Cachazo, Douglas, 
Seiberg, and Witten [6, 7, 8] pointed out that the planar diagrams describe the on-shell chiral 
ring of the gauge theory. They used a different technique to study the on-shell ring, namely 
the generalized Konishi anomaly equations, and those equations turned out to be exactly 
similar to the loop equations of the random matrix model. This confirmed the Dijkgraaf- 
Vafa gauge-matrix duality and made it more precise: the matrix model is dual to a subring 
of the gauge theory's chiral ring comprising the gaugino condensates and the mesons; other 
operators of the gauge theory are invisible to the matrix model. In particular, the sphere- 
level free energy of the matrix model is dual to the effective prepotential of the gaugino 
condensates. 

Gauge-matrix duality works for all kinds of theories, including quiver theories with mul- 
tiple gauge groups and bifundamental matter [2]. Here is a brief summary of Dijkgraaf- Vafa 
rules for building the matrix model of a particular gauge theory: 

1. A U{n) or SU{n) gauge symmetry of the field theory becomes a U{N) global symmetry 
of the matrix model, and we take the N —>■ oo limit. If multiple gauge symmetries 
SU{ni) are involved, we take all the iVj — oo at the same rate: Ni = t x rii for the 
same t — > cx). Similar rules apply to SO(n) or Spin) gauge group factors. 

2. Chiral superfields becomes bosonic variables of the matrix model in similar multiplets of 
the symmetry group. Thus quark and antiquark superfields become complex vectors 
of length iV, adjoint superfields become N x N matrices, and the bi-fundamental 
multiplets of an SU{n) x SU (m) symmetry become N x M matrices. 
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3. A complex chiral superfield becomes a real bosonic variable of the matrix model (or 
rather t or real variables. Thus a complex adjoint superfield becomes a hermitian 
N X N matrix. Likewise, a complex vector A of the matrix model and its conjugate 
B = together represent both the quark and the antiquark superfields of the gauge 
theory. 

Rule 3 works well for non-chiral gauge theories with real matter multiples (or conjugate 
pairs like □ + □) but not for chiral theories. For example, consider an SU{n) x SU (m) 
theory with a chiral bifundamental field (n, m) without the conjugate (n, m) multiplets.^^^ 
By rules 1 and 2 this field becomes an N x M matrix transforming as a bifundamental 
of the U (N) X U (M) symmetry, and since this representation is complex, the matrix must 
be complex too. But according to rule 3 such complex N x M matrix corresponds to the 
non-chiral (n, m)-|-(n, m) multiplet of fields rather than the chiral (n, m) without the (n, m). 

Lazaroiu et al. [9, 10] found a way out of this problem: instead of real bosonic variables 
integrated over R, one may use holomorphic variables, i. e. complex variables integrated over 
some contours in the complex plane rather than the whole plane C. More generally, one 
integrates N x M matrix elements of a complex matrix over some variety V C C^^ of real 
dimension NM. The variety F must be consistent with the symmetries of the matrix model 
according to rules 1 and 2. For non-chiral models one may identify F with the real "axis" 
■^NM (^NM recover rule 3, but chiral models require non-linear integration varieties. 
Lazaroiu et al. found such non-linear F for the model of a single-5'C/ {n) gauge theory with 
chiral B+(n — 4)0 matter spectrum [10], but in this article we interested in a different 
theory. 

Specifically, we are interested in the chiral bifundamental fields, hence we are looking 
for varieties oi N x M complex matrices which are invariant under the U{N)l x U{M)ji 
symmetry action. To be precise, the variety F should satisfy 

yVL^ U{N), yVne U{M) : V^TV^ ^ F (1.1) 

where '=' means equivalence as an integration variety: same topology with respect to the 

(1) The (n,m) bifundamental must be accompanied by other chiral multiplets with opposite anomaly, but the 
present argument does not depend on those multiplets. 
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singularities of the integrand, and similar asymptotics when one or more matrix elements 
approach the complex infinity. For square matrices N x N there is a simple solution, namely 
the unitary group space F = U{N) — which is actually invariant under the symmetry (1.1) 
(V^FV^ = F) and has the right real dimension = iV^. The problem is much more difficult 
for the rectangular matrices with M ^ N, and so wc leave them for future research. In this 
article, we work with the square matrices integrated over F — U{N). 

On the gauge theory side, a unitary matrix U E U (iV) is dual to an (n, n) bifundamcntal 
field Q with non-zero eigenvalues. That is, the vacuum states of the field theory which have 
matrix-model duals must have (det fl) ^ 0; there may also be vacua with (det O) = 0, but 
the matrix model might not work for such vacua. Indeed, under the gauge-matrix duality, 
the allowed values of the fields correspond to matrix variables belonging to the integration 
variety F or any of its allowed deformations F' = F. For example, a bosonic variable z 
integrated over a unit circle with measure ^ ^ can also be integrated over any other loop 
surrounding the z = origin. However, no such loop can go through the z = point itself 
because of the measure singularity, and consequently the field dual to z can take any 
complex values except = 0. Likewise, the unitary matrix integral can be deformed to 
integral over a variety F' = [/ (N) using an analytic extension 

Jdu;[U] = JUiU-'dU),, (1.2) 
r' r' (''^) 

pairs 

of the Haar measure. However, this extension becomes singular for det U — 0, which limits 
the F' integration varieties to the invertible matrices only. Hence, on the gauge theory side 
of the duality we should have det 7^ 0, otherwise the duality might break down. 

The bifundamcntal fields with detJl 7^ are common in dimensional deconstruction. 
Accordingly, in this article we build a unitary matrix model of an [SU{Nc)]^ 4D supersym- 
metric gauge theory which deconstructs the 5D SQCD [11]; the chiral ring of this theory was 
studied in great detail in [12]. The matter fields of the theory are shown on the following 
quiver diagram: 
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In particular, the blue lines here denote chiral bifundamentals, which we shall model via 
unitary matrices Ue, i = 1,2, . . . , K. The red lines denote the quark and the antiquark 
fields; we shall model them in the usual way as complex N x F matrices Ai and = A^. 

The article is organized as follows: in the next section §2 we build the matrix model 
of the gauge theory (1.3) and establish the gauge-matrix correspondence. In particular, we 
derive the the loop equations of the matrix model and see that they agree with the Konishi 
anomaly equations of the gauge theory. In §3 we evaluate the matrix integral and derive the 
free energy of the matrix model in terms of contour integrals on the spectral curve. We begin 
by reducing the problem to an integral over a single unitary matrix, and then we adapt the 
technology of Dijkgraaf and Vafa [5, 2, 1] and Cachazo et al. [6] to the unitary case. In §4 
we derive the effective superpotential for the off-shell gaugino condensates of the dual gauge 
theory. Finally, in §5 we discuss open questions related to the present research. 
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2. The Unitary Matrix Model 
and its Loop Equations 



2.1 The [SU{Nc)r Gauge Theory and its Matrix 
Model 

Dimensional deconstruction of 5D SQCD leads to an A/" = 1 [SU{Nc)]^ gauge theory in 4D 
as described in [13, 11]. The chiral ring of this [SU{Nc)]^ theory was analyzed in much 
detail in [12], and in this article we build and study its matrix model. We begin with the 
basic structure of the 4D field theory as shown in the quiver diagram (1.3): the green circles 
denote simple factors of the net gauge group 

K 

G4D = U^SUiN,)], (2.1) 

while the red and blue arrows denote the chiral superfields: 

^1 quarks J = (□<>), 

antiquarks Q{ = (D^), (2-2) 
*- link fields = (□f.i.i,^^), 

where / = 1,2, Nf and i = 1,2, . . . , K is understood modulo KS'^'^ Note that the link 
fields form chiral bifundamental multiplets of the gauge group. 

The tree level superpotential has three types of terms serving different purposes, 

Wtree = WqR + TThop + TTdef • (2.3) 

The O'Raifeartaigh terms 

Won = /^J^s, X (det(fi,) - i;^^) (2.4) 
i 

(2) For deconstruction purposes one takes the if — > oo limit, but from the 4D point of view K is a. fixed 
parameter of the field theory; in our analysis we shall assume K to be largish but finite. 
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— where the se are singlet fields not shown in the quiver diagram — turn each bifundamental 
field into an SL{Nc, C) linear sigma model. This is important for deconstruction purposes, 
and also allows us to model the fli with unitary matrices without worrying about the zero 
eigenvalues. The hopping superpotential 

w-hop = 7EE(^m^^^^./ - ^'fQ^Qi,f) (2-5) 

e=i f=i 

describes quarks' masses and interactions which let them 'hop' between quiver nodes; in 5D 
terms, this allows quark propagation in the deconstructed direction. Finally, we have the 
deformation superpotential 

d 

PFdef = tr(w{QK^K-i---Wi)) = ^^tr((f]i,---f]if) (2.6) 

p=i ^ 

for some polynomial W{X) — Ylp^^^ with constant coefficients gi,...,gd. The Wdef 
deforms the 4D theory away from deconstructed SQCD5, breaks the Coulomb branch of the 
moduli space into a discrete set of vacua, and allows formation of the gaugino condensates. 
This deformation is analogous to the tree-level superpotential for the adjoint field in [6, 14, 7] 
and is essential for understanding the on-shell chiral ring of the theory. It is also a key 
ingredient of the matrix model. 

In the matrix model, the [SU{Nc)]^ gauge symmetry of the field theory becomes a 
[U{N)]^ global symmetry and we take the TV — > 00 limit. Consequently, the bifundamental 
fields Qi become N x N unitary matrices Ue, and for chirahty's sake we should integrate 
each matrix Uc over the U{N) group manifold or an equivalent variety. However, in light of 
the O'Raifeartaigh terms (2.4) we restrict the determinants of the Ue matrices and integrate 
them over the SU (N) group manifold. To accommodate the v^" factor, wc rcscalc the field- 
matrix correspondence according to ^ v x where v = v + quantum corrections. Such 
quantum corrections are computable in the field theory — cf. §4.3 of [12] for details — but 
in the matrix model they need to be put in by hand. 

The quark sector of the field theory is non-chiral — for each quark there is an anti- 
quark Qi with opposite quantum numbers. In the matrix model, the quarks become complex 
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rectangular N xF matrices Ae while the antiquark fields Qi become conjugate F xN matri- 
ces B^ = Al. When we take the N ^ oo limit, we have two options for the fiavor number F 
of the matrix model: we may keep it fixed (i.e., F = Nf), or we may let it grow while 
keeping the fiavor/color ratio fixed, F/N = Nf/Nc [15]. In the 't Hooft limit of fixed F, 
the bifundamcntal sector dominates the matrix model in the large color limit, and the quark 
sector becomes quenched — its backreaction on the bifundamcntal sector becomes negligi- 
ble. This limit oversimplifies the physics but makes for a simple 1/N perturbation theory 
in terms of Feynman-like diagrams' topology. On the other hand, in the un-quenched limit 
of F, TV — > oo the matrix model has rich fiavor physics, but the 1/N expansion becomes 
much more difficult. Consequently, we use the 't Hooft limit in this article and leave the 
un-quenched fiavor physics for future research. 

Together, the Ui, A^, and matrices comprise the entire matrix model. Its partition 
function is defined as the following matrix integral: 

^ ^ / ^^[^^] J J dBe dAel cxp ( -jW {all Ue, Be, AeU (2.7) 

[sU{N) B,=A\ J ^ ' 

where the matrix potential is 

W{U,B,A) =W^,,{U) + W^o^{U,B,A) 



p=i ^ 

K 

-fJ2(vtT{Be+iUiAe) - tr(M^S^)) 



p 

K 

+ 

e=i 



where jj is an F x F matrix with eigenvalues fii, . . . , fip. Note that the matrix potential does 
not contain O'Raifeartaigh terms analogous to the FFqr of the field theory, but integrating 
the Uj( over the SU(/V") instead of U(/V") has the same effect. 

The denominator S in the exponent in eq. (2.7) is the overall coupling constant of the 
matrix model; as usual, S is fixed while N —>■ oo. Under gauge-matrix duality, S is dual to the 



8 



net gaugino condensate S — ■ Si of all subgroups of the SU (A^c)diag — diag[all SU {Nc)i] . In 
the field theory, the Si and hence the -S" emerge from the on-shell chiral ring, and then need 
to be integrated in to an effective off-shell supcipotential, but in the matrix model the S is 
an input parameter and the effective superpotential emerges from a more direct calculation 
we shall perform in §3. 

Finally, let us fix the overall normahzation factor C in eq. (2.7). Each Ue is integrated 
over a compact manifold, while each = Bj is integrated over the non-compact C^^, but 
the integral is Gaussian. Hence, we let 

C = X — 2.9 

(vol[SU(iV)])" 

where 1] is the Gaussian integral for a single quark mode of an average mass. The 'average' 
here does not have to be the arithmetic or the geometric mean, any representative value will 
do, and so we use ■jv because the modes vary in mass from ^{v — fif) to 7('0 -|- fif). Hence, 

r? = ldxdxexT)\ — —''yv\x\ = — — . (2.10) 



2.2 Loop Equations of the Matrix Model 

Having defined our matrix model we now need to verify that it is indeed dual to the 4D 
gauge theory which deconstructs the SQGD5. In this section we shall verify that the loop 
equations of the matrix model are similar to the anomaly equations of the gauge theory. To 
be precise, the matrix model is dual to a rather small part of the gauge theory — namely the 
subring of its chiral ring involving either the gaugino condensates or the mesons — but 
the anomaly equations for that part of the gauge theory should be accurately reproduced by 
the loop equations of matrix model. 

(3) The matrix models also has analogues of the baryonic and antibaryonic operators of the gauge theory, but 
such "dual baryons" are not matrices and are rather difficult to handle. For this reason, we shall limit our 
analysis here to the duals of the gaugino condensates and mesons only. 



9 



The loop equations follow from infinitesimal holomorphic changes of integration variables. 
For a simplified example, consider a toy model of a single unitary N x N matrix U, 



Z 



N 

duj[U] exp ( - J tT{W{U) 



(2.11) 



U{N) 



where (ia;[f/] is the holomorphic form (1.2) of the Haar measure. Let us change 



U ^ U' = U + ex f{U) 



(2.12) 



where / is a holomorphic function of U . Generally, this breaks the unitarity of U , hence we 
should deform the integration variety from T = U (N) to T' which spans the U' for U eU (iV). 
However, F' = F and hence this deformation does not affect the holomorphic integral (2.11). 
On the other hand, the variable change (2.12) itself has a non-trivial Jacobian 



J 



duj[U'] 
duj[U] 



{U-^dU\, 



jk 



dUi 



kj 



Nti[U-^f{U) 



(2.13) 



Also, the integrand of the matrix integral changes according to 



exp ||-^tr (!?([/)) j 



same x 



ejiv[W'{U)f{U) 



(2.14) 



Altogether, we have changed the matrix integral by 



5Z = eZx{{J2^-M^ - Ntr{u-^f iU)) - | tr(^-(^)/(^j )^ 
= e J duj[U] exp I ti (w{U)) ] x ^Vl 

U{N) 



(2.15) 



On the other hand, we have done nothing but changed the integration variable from U to U', 
hence the integral should not change at all. Therefore, for any holomorphic matrix matrix 
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function f{U) we must have 

(y^^-Mll _ Nir(u-^m) - jt^{w\U)f{U))\ = 0. (2.16) 

This toy example shows how to derive loop equations for unitary matrix models. Let us 
apply this technology to the unitary C/^ matrices of our big matrix integral (2.7). Let us pick 
one matrix, say C/^ and change 

^ U[ ^ + ex f{Ue, . . .) (2.17) 

where / is a holomorphic function of the Ui and other matrices of the model (denoted 
by the . . .). To preserve the symmetries of the model, f{Ue,...) must be covariant, i.e. 
transform like Ui under the [U{N)]^. Also, because the is integrated over the SU(iV) 
group manifold rather than U(7V), we must preserve the det(C/^') = det{Ui) — 1 condition, 
hence / must satisfy 

tr (u^^fiUi,...)) = 0. (2.18) 

On the other hand, we do not need to preserve the unitarity of the U'^ because we may 
deform the integration variety from F = SU{N) to a nearby F' = F. Hence, proceeding 
exactly as in the toy example above, we find that for any covariant holomorphic / which 
satisfies eq. (2.18) we must have 

Now let us focus on functions / which depend only on the unitary link matrices but not 
on the quark matrices. By covariancc, products of the Uei must be taken in the order of the 
quiver, hence f{Ui, other ^7^/) must be a hnear combination of (t/^t/^-i • • • UiUk ■ ■ ■ Ui+iyU^ 
for p = 0, 1, 2, . . ., which can be summarized in a power series 

(2.20) 

in an auxiliary complex variable X. In field theory we use this series as it is, but in the 
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matrix model we must correct for the trace condition (2.18), thus 



f{Ue,...) 



X - f/,f/,_i ■ ■ ■ UiUk ■ ■ ■ U,+i S 



X Ui (2.21) 



where the second term inside the brackets assures tr[*] = 0. Specifically, 



^^^^ N^'^^X - U,U,.i ■ ■ ■ UiUk ■ ■ ■ U,+i 



(2.22) 



which is same for all i and remains finite in the N ^ oo limit. For / as in eq. (2.21) 



E 

jk 



dfkjiUi 



XR\X) - SR{X) + -^{XR'{X) + R{X)) 



(2.23) 



while 



where 



XW{X)R{X) - SW{X) + SP{X) - CR{X) 



(2.24) 



W'(X) 
P{X) 

and C 



p=i 

— tr [Uk ■ ■ ■ Ui) 

N V X - Uk---Ui 

(a polynomial of X of degree = d — 2) 
J ti({UK---U^)W{UK---U^)). 



(2.25) 
(2.26) 

(2.27) 



Note that in the large limit C and P{X) remain finite, hence the right hand side of 
eq. (2.24) grows like A^. By comparison. 



trl^m,...)) =0(F), 



(2.28) 
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which is much smaller than 0{N) in the 't Hooft limit. Therefore, eq. (2.19) becomes 
'XR\X) - {XW'{X)+S-C)xR{X) + SW{X) - SP{X)) = 0{1/N) ^ 0. 



(2.29) 

Moreover, for N ^ oo matrix averages factorize as 

{R\X)) ^ {R{X)f , {C X R{X)) ^ {C) X {R{X)) , (2.30) 
and this gives us a loop equation for the {R{X)), namely 
X{R{X)f - (^XW'{X) + S-{C)^ X (i?(X)) + SW'{X) + S {P{X)) = 0. (2.31) 

On the field theory side, we have a similar quadratic equation for the gaugino condensate 
resolvent 

'^f^') = ^ .Y-n,n._...'.n.n....n,,. ) (^-'' 

(2.32) 

In the on-shell chiral ring, this resolvent satisfies 

XR\X) - W{X)R{X) + F{X) = (2.33) 

where 

d 

VV(X) = J^dp^^ + P^^^(^) (same {sg) V£) (2.34) 

p=i 

and F{X) is another polynomial (of degree < d) which depends on vacuum state of the field 
theory. Eqs. (2.31) and (2.33) are obviously dual to each other; to make them identical we 
simply need to identify 

X^ly — X X-tr, RiX)'lZry " ^^"^ X (/^(X))-^^!^ , (2-35) 

iPv'''{s))l^eU ^{S- (C))rder, HX)^Zry — ^)-^ X ^ (W'(X) + {PiX)))^^^ . 

(2.36) 
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This correspondence explains why S is dual to the net gaugino condensate 



. . (2,37) 



Solving the anomaly equation (2.33), we have 



^ W(.Y) - VW^(.Y) - 4A-F(A-) _ ^ ^ 
2A' 

According to eq. (2.36), the matrix dual of the right hand side is 

AH"(A)+ const X^Sf A'S,, ^ ^ 

(note that P{X) has lower degree in X than W(X)), and therefore 



('S'net) theory ^ ('S')modcT- (2.40) 

The mesonic resolvents of the gauge theory also have matrix duals. On the field theory 
side, chiral mesonic operators with quarks and antiquarks at different quiver nodes are 
packaged into a bunch of resolvents 

Me'AX) = Qi'^e'-i ■■■^e ^ ^ p- Qi (2.41) 

J\ — • • • ■ ■ ■ iZ^ 

subject to periodicity conditions 

-Mi>=t+K + XMi>=t = Me = {QiQt) (2.42) 

where the right hand side is an ordinary mesonic operator which does not depend on X. 
Besides their quiver indices, the mesonic resolvents are also Nj x Nf matrices in the flavor 
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space. On the matrix-model side, they are dual to (the averages of) F x F matrices 

M,,e(X) = BeUe-i---Ue- ^— -A, (2.43) 

A — Ui-i ■ ■ ■ UiUk ■ ■ ■ 

which satisfy a similar periodicity relation 

-Me,=i+K(X) + XMe=e{X) = M, = (S,A,). (2.44) 

Because of the v factors in the <-> C/^ x correspondence, in the mesonic sector we expect 
the gauge-matrix duality to work according to 

^tSy - X Xr^, M,AX)t^ry - v''-'-'' X {M,,e{X))-^^ , (2.45) 

To derive loop equations for the mesonic resolvents of the matrix model, we need to 
vary the and the matrices independently of each other. Note that the conjugacy 
constraint = is just a special case of integrating the 2NF complex numbers comprising 
each {Ai, B() over a variety V of real dimension 2NF, we "just happened" to choose F = 
{B^ = ^4^}. But as long as the integrand is a holomorphic function of both Ai and B£ 
(views as independent variables), we may deform the integration variety without changing 
the integral, 

JJ dAe dBefiAe, Be,...) = JJ dA, dB, f{A,, B,, . . .) for V ^ {B, = A\}. (2.46) 

Therefore, sma// variations of the Ap^ and the B^ matrices don't need to be conjugate to each 
other — the discrepancy will deform the integration variety a bit, but small deformations 
do not affect the integral. 

Minding this rule, let us vary any one B^ matrix while the "conjugate" A(, matrix remains 
unchanged. Specifically, let 

B, ^ B[ ^ B,+ exBt.x ^^''-\"'^\^ (2.47) 
where e is an infinitesimal F x F matrix in the flavor space. The Jacobian for such change 
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of variables is 



1 + Se',extr{e)xti(- -i -) = 1 + Se,extr{e)x^ R{X), (2.48) 



while the potential changes by 



5W = ^T^hop = IV ti (e X Me,e-i{X)) - 7tr (e x M,,,,(X) x fi) . (2.49) 

But altogether, this is just a change of an integration variable which does not change the 
integral at all, hence 

= (5(Jacobian)) - ^ im\ (2.50) 

S ^ ' 

= X tr(e) X ^ (i?(X)) - ^^tr (e x (M,,,,_i(X)) - M,,,,(X) x /i)) . 
This must hold true for any e matrix, therefore 

V (M,,,,_i(X)) - (M,,,,(X)) X /i = V X 7"' (i^W) X 1^^^ . (2.51) 

Likewise, changing the An matrix by 

?>A, = X X e (2.52) 

while the Bi matrix remains unchanged leads to another loop equation, namely 

V (M,+i,,,(X)) - /i X (M,,,,(X)) = V X (^(^)) X 1^^^ . (2.53) 

Both loop equations (2.51) and (2.53) look exactly hke their field theory counterparts: 

Mi>,i-x{X) - Me',e{^)xfjL = Me'+i,e{^) - pLxMe',e{^) = V x 7-^ i?(X) x Ijv^xJVy 

(2.54) 

The only difference is in the v factors — and that is in accordance with the correspondence 
rules (2.45) and (2.35). 
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Having two sets (2.51) and (2.53) of mesonic loop equations is not redundant. Combining 
both sets with the periodicity equations (2.44) allows us to solve for all the mesonic resolvents 
(or rather their averages) in terms of a single matrix (M) . The solution works exactly as for 
the gauge theory, so instead of copying from [12] almost verbatim, let us simply state the 
result: 



(M,,=,(^)) 



(M,,>,(X)) 



(M) + {R{X)) X 1 



'-^-1 /i X (M) + {R{X)) X 1 



X 



X - v-Kfi^ 



(2.55) 
(2.56) 



where (M) = {BpAp) must be same for all ^ and must commute with the quark mass 
matrix /i.^^^ 

Finally, there one yet another bunch of loop equations for the mesonic resolvents stem- 
ming from A and B dependent variations of the U matrices. Let us pick any one Ui and 
vary it according to 



5Uf 



[Ai+i X ex BiX 




X - U,.i---UiUk---U, 



X Ui (2.57) 



where the second term assures 5(lei{Ui) = 0. The trace in the second term evaluates to 



tr e X 



M 



e+N,e+i 



(X) - Me+N,e+i{0) 
X 



{2.5t 



The Jacobian of the variable change (2.57) is given by 



5j = j: 



d{5U, 



e)jk 



jk 



e)jk 



N 

— tr 

S 



e X 



V 



■R{X) X Me+K,i+i{X) 

Me+K/+i{X) - Me+K,e+i{0) 



S X 



L + 0(l/iV2) 



X 



(2.59) 



(4) Strictly speaking, (M) must commute with the fi^ matrix rather than with the /t itself. Also, the matrix 
elements of {B(Ai) which commute with fi^ but not with jl have ^-dependent phases. To avoid this mess, 
we assume that all distinct mass eigenvalues also have distinct K— powers. 
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while the potential varies according to 



5W^hop = 7«tr(exM,,,(X)xM,+i)+ 0(l/iV), 



SWdef = = tr 



e X 



v 



■W'{X) X Me+KMX) 



- C X 



tiiX) 



X 



J/ 



(2.60) 
(2.61) 



where C is as in eq. (2.27) and t£{X) is an F x F matrix-valued polynomial of X whose 
actual form does not affect the following argument. Altogether, these changes must cancel 
out of the matrix integral, hence demanding 



N 

{5J) - -^{SW) = Ve 



(2.62) 



we arrive at 

{Mt,e(X) X M,+i) - {R(X) x Me+K,e+i(X)) 

S-C 



W{X) + 



X 



S-C 
X 



UX)) + 0{1/N). 



(2.63) 



In the large N limit matrix averages factorize, and this gives us another family of loop 
equations, namely 



jv{Me,e{X)) x (M^+i) + 



W{X) + 



S-C 
X 



{R{X)) 



X {Me+K,e+iiX)) 



S-C 
X 



(2.64) 



a polynomial of X. 



Again, these loop equations are exactly similar to their field theory counterparts and we 
may solve them in a similar way. Making use of eqs. (2.55-56) and eq. (2.31), we transform 
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eq. (2.64) into a quadratic equation for the mesonic "VEV" (M), namely 

(li^mf + |^l^W'(X) + S-Clx(T/i(M)) + ^F(X) 
~ \ ^ ) ^ polynomial of X) 

Also, the (M) matrix is block-diagonal in the eigenbasis of /t and does not depend on X. 
Hence, for each block we may substitute a different value of X into eq. (2.65) and apply 
the resulting equation to the block in question; our choice \s X = jlj /v^ which kills the 
right hand side of eq. (2.65) regardless of the matrix-valued polynomial we didn't spell out. 
Consequently, each block — and hence the whole matrix — satisfies 

(7/i(M))2 + ((/i/«)^>V'((/i/i;)^) +S-C)x (7/i(M)) + F((/i/t;)^0 = 0. (2.66) 
Note similarity between this equation and the loop equation (2.31); this allows us to write 
(M^) = -^MAIL evaluated for X = (/i/i))^ (2.67) 

where {R{X))^ on the right hand side indicates the two solutions of eq. (2.31). 

This completes our analysis of the loop equations of our matrix model. Having seen 
that those equations are dual to the anomaly equations of the gauge theory of the quiver 
diagram (1.3), we can be positive that our model is indeed dual to that gauge theory. 
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3. Calculating the Matrix Integral 



In this section we evaluate the matrix integral (2.7) in a sequence of simple steps. First 
we integrate over the Ai and Bi matrices dual to the quarks and the antiquarks. Second, 
we reduce the integral over K link matrices Ue to an integral over a single unitary matrix 
U = UkUk-1 ■ ■ - Ui. After that we follow the Dijkgraaf-Vafa method adapted to a unitary 
rather than hermitian matrix: we reduce the integral over the whole matrix U to an integral 
over its eigenvalues, and then we use the saddle-point approximation in the large N limit. 
In this limit, the eigenvalue spectrum becomes continuous with density p(A), the free energy 
has a 1/iV expansion where the sphere-level and the disk-level terms are given by spectral 
integrals, — and we relate the whole shmeer to the loop equation (2.31) and the period 
integrals of its Riemann surface. 

Let us integrate over the quark matrices. The hopping potential (2.8) is bilinear with 
respect to the matrix elements of the and Bf 

= Y:>y£B,,,,v'%A^-f (3.1) 

hi' f,f' 

where V'''% = v6%, (Uf, 5^ - 5f 5^ /i^ . (3.2) 

This makes the quark-matrix integral a Gaussian integral over KNF independent complex 
variables, which we evaluate as 



K 



i=i 



)KNF 



B,=4 



where the determinant involves all indices: quiver, color, and flavor. In block form 



V 



( ~A 








vUk 


vUl 













VU2 


-fl 






(3.4) 



\ ••• vUk-i -fl ) 
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and hence 

F 



Det {V) = Det (v^UkUk-i ■■■Ui - fi^) = TT det [v^UkUk-i • • • C/i - /tf ) 

KNFxKNF NFxNF J-_^NxN 

(3.5) 

where on the right hand side the /iy are eigenvalues of the quark mass matrix fi. 

Note that the 'quark' integral (3.3) depends on the link matrices Ui only through their 

product U = UkUk-i ■ ■ - Ui. Likewise, the deformation superpotential Wdcf depends on the 
only through their product. Hence, at this stage, we may write the matrix integral (2.7) 

as 



n ] / ^^[^^] \ = UkUk-1 ■■■Ui) (3.6) 



(vol[SU(7V)]y 

where 

f{U) = exp ( -^WaAU) ] I n det [u - (////i))^), (3.7) 

and d'jj\^(\ is the Haar measure for Ui e SU{N). This measure is left-invariant — for any 
fixed V e SU{N), du[VU] — duj[U] — and this makes it easy to change variables in the 
unitary matrix integrals. In particular, in an integral over K matrices such as (3.6) we can 
set Ui — VU{ (for £ — 1 only) and have 

dupK] X • • • X du;[U2] X dupi = VU[] = dupK] x • • • x du;[U2\ x du;[U[] (3.8) 

for any V e SU{N) which does not depend on the U[ matrix, even if V depends on the 
other unitary matrices U2, ■ ■ ■ ,Un- Therefore 

J duj[UK] ■■■ J duj[U,]f{UK^^^U,) 

SU{N) SU{N) 

= j dupK] ■■■ j du[U2\ j du[U[\f{UK^^^U2^VU[) 

SU{N) SU{N) SU{N) 

((setting V = (C/^ • • • ^^2)"^ ^ y-^C/i ^ Uk ■ ■ ■ U2U1 = U)) 



j dujpK] ■ ■ ■ j du[U2] j du[U[ = U] f{U) 



SU{N) SU{N) SU{N) 
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= (vol[SU(iV)]) 



K-l 



X / duj[U]f{U), 

SU{N) 



(3.9) 



and hence 



Z 



Vol[SU(Ar)] 



X 



duj[U]fiU). 



(3.10) 



SU{N) 



Furthermore, according to eq. (3.7), f{U) depends only on the eigenvalues of the unitary 
matrix U but not on its eigenvectors. Consequently, we decompose the — 1 coordinates 
of the SU{N) group manifold into N — 1 independent eigenvalues e*'*'' {i = 1, . . . ,N but 

Xi = mod 27i) and N{N — 1) angular variables 9i^j describing the eigenvectors. The 
Jacobian of this decomposition is given by 



dfOO 1 — r o Xi — 

— — — — = const X 4 sm^ — — ^ 



(3.11) 



— the unitary version of the Vandermonde determinant — which depends only on the 
eigenvalues Aj, hence integrating over the eigenvector variables 9ij we arrive at 



-,KNF 



2n 



2n 



dXi f dX 



N\ J 2tt J 2tt 





/(Ai,...,A^) X Y[Asm 



2 Aj — Xj 



X 



i<j 



L=—oo 



(3.12) 



where the last factor is the 5-function for the J2i modulo 27r. 

Thus far we made only exact calculation, but now we turn to approximations valid in 
the iV — > oo limit. Re-writing the integrand of eq. (3.12) in exponential form 



exp 



-^^W(e^^O -EElog(e^^-(/.//^)") 

.7 = 1 



3=1 f=l 



i<j 



:Sm 



N 



(3.13) 



we see that all terms in exponent grow with N. Hence, in the large limit we may use the 
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saddle-point approximation: 



^ N N F 



logZ ^ -jJ2^{e^''^ -J2T.^^^{''''-(l'f/^f) (3-14) 

j=i i=i /=i 

+ 5;iog(4sin^^) + ilf^A, 

i<j ^ ^ j=i 

where (A^, . . . , A^) and L maximize the right hand side of this formula, or rather maximize 
its real part and extremize it imaginary part; generally, this requires moving away from the 
real axis into the complex plane. Or into some other complex space: since the Aj are periodic 
variables on the circle ]R/27rZ, the Aj move into the complex cylinder C/27rZ rather than 
into the plane C. As to the L, in the large N limit the maximum happens for L — 0{N) 
where the discreteness of L does not matter any more, hence L — ^ x C where C is finite 
and complex. 

As usual for matrix models, for iV ^ oo the spectrum E of (Ai, . . . , A^) becomes contin- 
uous. In general, it comprises several continuous line segments Ei, . . . , E„ on the complex 
cylinder. Denoting the spectral density over those lines by ^p(X)dX, we write the free energy 
of the matrix model as 

:F = -S- logZ = ^5 + i:Fn + 0(1/N'), (3.15) 
iV"^ J\ 

where 

^ J dXp{X) W{e'^) - iCX - \ j dXp{X) J dX' p{X!) \og{^sui^^^^ (3.16) 



whole whole whole 

S EE 



and 



^ r 



f =1 
•' whole 

E 



In "worldsheet" terms of the 1/ N expansion [14], the leading contribution which comes 
from Wdef and the Vandermonde determinant corresponds to genus g = i.e. spherical 
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topology, hence the notation. Likewise, the quarks' sector contribution corresponds to the 
disk topology, hence the notation Tb. 

The spectral density p(A) minimizes the free energy T ~ J-s under constraint 

#(A0 ^ J dXp{X) = S, (3.18) 



whole 
E 



and this gives us a variational equation 

d 5Ts 



d\ Sp{\) 

In light of eq. (3.16), this equation becomes 



0. (3.19) 



-dx - 

whole 
E 



J dAV(A') cot = (3.20) 



where the bar across the integral sign indicates principal- value integration over the pole at 
A' = A. However, it is more convenient to regulate the pole via bypassing it in the complex 
plane, so let us introduce the resolvent 

R{w) = /rfAp(A) cot^^^^ (3.21) 



whole 
E 

which has branch cuts along the spectral segments Ei, . . . , S„ but is analjd;ic elsewhere in 
the complex cylinder. In terms of this resolvent, eq. (3.20) applies to the average of the two 
sides of any branch cut, hence 

VAeS: R(X±ie) - iC ^ 2mpi\). (3.22) 

dX 

Together, eqs. (3.21-22) lead to a quadratic equation for the resolvent R{w): 

'dWie'"^) 



R\w) - 2 



dw 



-iC 



X R{w) + 2F(e'"') + 5^ = (3.23) 
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where F is a polynomial function of e*"'.^^^ Physically, eq. (3.23) is nothing but the loop 
equation (2.31) written in terms of the periodic coordinate w of the complex cylinder instead 
of the fiat coordinate X — 

{R{X)) - 

and hence 

R{w) 

Likewise 

and if we also identify C — (C) then eq. (3.23) becomes the loop equation (2.31). 
Solving the loop equations gives us 

{R{X)) = X^'{X) + S-C - Y ^ ^ ^ iXW'iX) - id - iY{X) (3.28) 

zX 

where 

Y{X) = ±^(XW{X) + S-Cy - 4SX(W'{X) - (P(X))) (3.29) 

has n < d — degree(W) branch cuts which connect simple zeroes of the discriminant Y'^{X). 
Let us denote such zeroes a~ and (i = 1, . . . , n)^^^ and map them into the w cylinder as 



Indeed, in the saddle point approximation we have 



N X - exp(iAi) J X -e<-^ 

* whole 

S 



2iX {R{X)) - iS for X = e^"'. (3.26) 



= ^ X ((C) - XW'{X) + 2X {P{X))) , (3.27) 



(5) Integrating / d\' p{\') J dX" p{X") the trigonometric identity 

\±ie-\' \'±ie-\" , \±ie-\" ( A±ie-A' A' ± ie - A'' 

cot X cot 1 = cot X cot h cot 

2 2 2 V 2 2 

we arrive after some work at eq. (3.23) where 

FH =ydAp(A) cot^ X [-^ - ^-^j . (3.24) 

E 

(6) From now on we use index i to label spectral segments Sj and related parameters rather than individual 
eigenvalues. 



25 



af = exp{i9^). Consequently, the A spectrum S has precisely n segments Sj which begin at 
9~ and end at 9^', the spectral density along each segment is given by 



P(A) = ^Y{X = e^') 



(3.30) 



(it cf. eqs. (3.22) and (3.28)), and the shapes of the segments follow from p{X)dX being real 
and positive. 

From the field theory point of view, the segments Sj correspond to confining (or rather 
pseudo-confining) subgroups SU{Ni) of the diagonal SU{Nc), and the integrals 



ef 

Si = ld\p{\) = ldwR{w) = ldX{R{X)) 
J 4:711 J zm J 

g- Ai Ai 



(3.31) 



are dual to the gaugino condensates 



S,. = Tr 



327r2 



SU{N,) 



— (pdXR(X). 
2m J 



(3.32) 



In both cases the are A-cycles of the Riemann surface of the -R(-^)fheory ^ (-^("^))modeT 
or their images on the w cylinder; the Ai surround individual segments of the spectrum as 
shown below: 





X = 





(3.33) 
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The advantage of the contour integral formahsm is that we do not need to know the exact 
routes of all the segments but only the general locations of the ^-cycles which surround them. 
This works for all integrals of the form JdXp{\)f{\), for all kinds of /(A) functions. For 
example, eq. (3.17) for the disk-level free energy can be written as a sum of contour integrals 



(3.34) 



The sphere-level free energy can also be written as a sum of contour integrals, but this is 
more complicated because the logarithm in the second terms in cq. (3.16) has branch cuts 
of its own. To disentangle the branch cut structure, let us first rewrite eq. (3.16) as 



^ \ j d\p{\) X W{e'^) - iCX + E{X) 

whole 



(3.35) 



where 



Naively, 



E{X) = W{e'^) - iCX 



- JdX'piX') 



X log ( 4 sin 



A- A' 



whole 
S 



MA) " ^ ^ ^ rfA - ° 



[cf. eq. (3.19)), but analytic continuation to the complex cylinder yields instead 



(3.36) 



(3.37) 



dE 

Ix 



dX 



iC 



R{X) = iY{X^e'^) ^ 



(3.38) 



and hence the E{X) itself is a non-trivial function with branch cuts along the Ej segments. 
However, its average between the two sides of a cut is locally constant 



^(A + ie) + E{X - ie) 



constant Ei for A e Ej only 



(3.39) 



because the derivative (3.38) flips sign across the cut. Therefore, to make sure eq. (3.19) is 
consistent with the loop equation, we must disambiguate the logarithm in the second term 
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of eq. (3.16) such that 



for A e E, 



5J^s E{X + ie) + ^(A - ie) 



5p{\) 



Hence in hght of eq. (3.39), 



(3.40) 



i=l 



JF^ = d\p{\) X \w{e'^) - iCX + E, 



1 " 



i=l 



EiXSi + — IdwRiw) X (wi^'^) - iCt 
Am J \ 



(3.41) 



The Ei can also be calculated as contour integrals. For general A, 

E{\) = W(e'^) - iCX - jdwR{w) x log Usm^^^^\ (3.42) 

where the contours should be drawn such as to exclude the branch cuts of the logarithm. 
In particular, the A point itself should be kept outside of all the contours, and this prevents 
us from directly evaluating eq. (3.42) for A e E. Instead, for the purpose of calculating an 
Ei in eq. (3.41) we must first evaluate the integral eq. (3.42) for A ^ E, then take two limits 
of A approaching the same point Aj e Ej from two opposite sides of the spectrum, and finally 
take the average of the two limits. Alternatively, we may take just one limit of A going to 
an end point Of or 6^ of the spectral segment — at these points the difference between the 
two sides of the spectrum vanishes and the averaging becomes unnecessary. Thus, 



E, 



lim 



W{e'^) - iCX 



1 

- —Y 



dw R{w) X log ( 4 sin 



X — w 



(3.43) 



In the next section, these equations (as well as eq. (3.34)) will help us calculate the effective 
superpotential of the matrix model. 

We conclude this section with a brief discussion of local versus global minima of the free 
energy. Eqs. (3.41) and (3.43) give us the J-'s minimized with respect to local variations of the 
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spectral density p(A), hence eq. (3.30). The global minimum requires further minimization 
with respect to the free parameters of eq. (3.29), namely C and n — 1 independent coefficients 
of the (P(X)) polynomial. ^'''^ Variationally, this implies that 5J^s/^p{^) should be globally 
constant over the whole spectrum S and not just individual segments, hence 

ii = 4 = ••• = (3.44) 
In addition, we should minimize with respect to the C parameter, hence 

5:Fs 



sc 



j d\p{\) X {-iX) = -^X] jdwR{w) xw = 0. (3.45) 



vhole Ai 
S 



Note that eqs. (3.44) apply only to the global minimum of the free energy; in field theory 
terms this corresponds to taking the gaugino condensates Si on-shell. In the following section, 
we shall calculate the Weff for the off-shell Si or rather S^, and this means abandoning the 
global minimum and hence eqs. (3.44). Instead, eqs. (3.31) will determine the coefficients of 
(P(X)) in terms of the Si. 



(7) {P{X)) has degree d — 2 and hence d—1 coefficients, but if we want n < d spectral segments, the polynomial 
Y'^{X) (cf. cq. (3.29)) must have d — n double zeroes, which imposes d — n constraints on the coefficients of 
{P{X)). Consequently, only n — 1 of those coefficients may vary independently of each other. 
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4. The Effective Superpotential 



In this section we derive the effective superpotential of our matrix model. But first, a few 
general words about effective superpotentials for the off-shell gaugino condensates in field 
theory. In the single-f/(A^) theory with adjoint matter, Cachazo et al. [6] found that the 
gauginos of the U{1) center of the U{N) generate auxiliary supersymmetries of the chiral 
ring of the gauge theory and hence 

1=1 ' i=l 

where J-'{Si, . . . , S'„) is a prepotential of the auxiliary SUSY, tq is the overall bare gauge 
coupling, and bi are integers distinguishing between specific vacua of the theory. Under 
gauge-matrix duality, tq is just an arbitrary parameter, but the prepotential J-' is dual to 
the sphere-level free energy J^s{Si, ■ ■ ■ , Sn) of the matrix model. 

Adding the quarks — which couple to the U{1) center — breaks the auxiliary super- 
symmetries, but the subring of the adjoint sector remains super symmetric. Hence [14], 
integrating out the mesonic sector of the chiral ring yields 

n n 

W,s{Si,...,Sn) = J]iV,— ^(^i,...,^^ +J2i^^^^o + k)xS, + J^Q{Su...,Sn) (4.2) 

i=l * i=l 

where J-'q is the quark sector's contribution to the superpotential. Under gauge-matrix 
duality, the J^q is dual to the disk-level part J-'d of the matrix model' free energy. 

In our case, the [SU{N)]^ theory does not have a U{1) center(s), and we cannot promote 
the gauge symmetry to [U{N)]^ because of gauge anomaly constraints, but we can add a 
single U{1) factor common to all nodes of the quiver diagram (1.3). Since all the gaugino 
condensates of our theory belong to subgroups of the diagonal SU{N), adding the common 
f/(l) effectively promotes the SU{N)diag to the f/(A^)diag, with similar consequences for the 
effective superpotential. Thus, Wes{Si, . . . , Sn) is given by eq. (4.2) where the J-'{Si, . . . , Sn) 
is dual to the J^s of our matrix model and the J-'q{Si, . . . , Sn) is dual to its J-'d- 

(8) By analogy, in the Veneziano-Yankielowicz superpotential W — NSlogS + {2mTQ + 6)5', h distinguishes 
between different vacua of the SU{N). 
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In the matrix model, taking the Si off-shell and evaluating the free energy as a function 
^^(5*1, . . . , Sn) means imposing all n eqs. (3.31) as constraints and then minimizing the ~ 
J-'s under those constraints. Variationally, this implies that 6J-'s/Sp{X) should be locally 
constant along each continuous segment Sj of the spectrum but may take different values for 
different segments. Therefore, eq. (3.30) and all the subsequent contour-integral formula of 
§3 remain unchanged, but we throw eqs. (3.44) out of the window and instead use eqs. (3.31) 
to determine the coefficients of (P(X)). Consequently, exactly as in [6] 



= Ei, (4.3) 



and hence the effective superpotential of the matrix model is given by 

n n 

W,s{Si,...,Sr.) = 5^iV,xi,(^i,...,^„) +5^(27rzro + 6.)xS, + J^d{Si, . . . , S^)- (4.4) 

i=l i=l 



In the contour integral formalism, the Ei in this formula are given by eqs. (3.43). We 
can simplify the contour integrals by merging all the Ai cycles together and then pushing 
the resulting loop outwards. The following picture illustrates how this works on the complex 
cylinder: 




The new contour comprises two vertical lines on two sides of the logarithm's branch cut, and 
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two loops around the cylinder. For the loops, we take Q — > oo which gives us 



R{w) 



Im(u>) ±00 



and therefore 



(4.6) 



top loop = 27riS (Q — Im(Ai)) , 
bottom loop = 27riS (Q + Im(Aj)) . 



For the vertical lines, we have discontinuity 

\ — w 



disc 



log I 4 sin 



+2-111 for Imw > ImA, 
— 27ri Imw < ImA, 



and hence 



dwR{w) X log(- • •) 
Putting the whole contour integral together, we arrive at 



vert, 
lines 



2ni J dw R{w) + 2m JdwR{w). 

A A 



+iQ —in 

^(A) = W{e'^) - iCX - ^ ldwR{w) - ^ ldwR{w) - SxQ 



A 

+'too 



W{e'^) - iCX - - I dw[R{w)+iS] - - I dw[R{w)-iS], 



(4.7) 
(4.8) 



(4.9) 



(4.10) 



(4.11) 



and at this point, we may take A = 6*^ or A = 6^ to calculate the Ei. 

We may also express the in term of the 5-cycle periods of the Riemann surface of 
Y(X — e*"'). Indeed, plugging eq. (3.28) into formulae (4.11), we obtain 



Ei ^ ^- J dwYie'"^) + ^- J dwYie'"^) + iW(e*"') - Sfl 



rcg 



reg 



dwYie'"") - - I dwY{e"") + i>V(e™) - SQ, 



(4.12) 
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where the Bf cycle begins at w = ±ioo on the physical sheet of the Riemann surface, crosses 
the Sj branch cut to the other sheet, and then goes back to w = ±ioo but on the unphysical 
sheet, and the integrals are regularized by starting and stopping at w = ±ifl instead of 
w = ±zoo. The cycles are illustrated on the following figure: 



Imw = +Q 




(4.13) 



Imw = —Q 



Note that the each sheet of our Riemann surface is a cylinder with two distinct infinities — 
and that's why we have a double set of -B-cycles. Fortunately, most of the extra B cycles 
are redundant: 



: B^ — B- = Bf — B- modulo A cycles. 



(4.14) 



and hence 



reg 



-- IdwYie'' 



+ terms common to all Ej 

+ a linear combination of Sj with integer coefficients. 



(4.15) 



The last term here reflects the ambiguity of 5-cycles modulo A cycles: going from a branch 
cut Sj to ±zoo one may choose different passages between the other branch cuts S^. In 
terms of the sphere-level free energy, this corresponds to different routing of the branch cuts 
of the log (4 sin^ around the spectral segments Sj in both A and A' complex cylinders. 
Re-routing the log's branch cuts changes the free energy by 



(4.16) 
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for some integer coefficients Qj, and hence 

E,{Su...,Sn) ^ MSu...,Sn) +J2c^,S,. (4.17) 

j 

Fortunately, from the superpotential (4.4) point of view, this change amounts to changing 
the integers foj. In other words, we permute the branches belonging to different vacua, but 
the overall picture does not change. 

Now consider the quark sector's contribution J^o{Si, ■ ■ ■ , Sn) to the effective superpoten- 
tial. Again, we are going to simplify eq. (3.34) by moving the integration contours, but it's 
convenient to do it flavor- by-flavor. We need to distinguish between massive and massless 
flavors: they play different roles in fleld theory, both in 4D and in 5D^^\ and in the matrix 
model they have different contributions to the disk-level free energy. Let us begin with a 
massive flavor and let e^'^f = {jlf/v)^. Then for this flavor we have cuts and contours as 
shown on the left picture below: 




(4.18) 

Again, we merge the Ai contours into a single loop and push it outward as shown on the right 



(9) In deconstruction [11], massive 4D flavors with k, v have Hght modes and give rise to hght flavors in 5D; 
the massless 4D flavors do not have no 5D counterparts, but they are needed to adjust the Chern-Simons 
level of the 5D theory. From the purely 4D point of view [12], the massless flavors affect the Coulomb branch 
of the theory but do not give rise to Higgs branches; the on-shell mesonic and baryonic operators include 
only the massive flavors. 
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picture above. This time, we end up with two disconnected contours: a simple loop around 
the top of the cylinder, and a loop around the cylinder's bottom attached to a vertical loop 
around the logarithm's cut. The discontinuity across this cut is simply +2m and hence 



^dwR{w) X log {e'"" - e'^^f) = 2m j dw R{w). 

-in 



(4.19) 



vert, 
loop 



And for the loops around the cyhnder we again take f2 ^ oo and use eq. (4.6); consequently 



top loop = 2mS X {imf — ni), bottom loop = 2TTiS x {Q + ni). 



(4.20) 



Putting all the loops together, we find that one massive quark flavor contributes 



I f . . S . . 1 



dw [R{w)+iS] +iSxmf. (4.21) 



-in 



For a massless flavor we have log (e™ — (O/'O)'^) = iw, which does not have singular- 
ities on the complex cylinder but needs a branch cut anyway because w is multi-valued. 
Consequently, we have contours and cuts as shown below: 




\m.w = —VL 



m • 






(4.22) 

Once again we merge the contours into a single loop and push it outward as shown on the 
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right picture above. This time we get a connected system of two loops around the cyhnder 
and two vertical lines, and evaluating the integrals gives us 



vertical lines : =^ 27ri JdwR{w), (4.23) 

-in 

each horizontal loop : =^ 27TiS x {fl ± ni), (4.24) 



hence altogether one massless flavor yields 

+ioo 



^D\Z:r" = I Jdw[Riw) + iS] + \ jdw[R{w)-iS]. (4.25) 



-lOO 



We conclude this section with a complete formula for the effective superpotential. Com- 
bining eqs. (4.4), (4.11), (4.21), and (4.25) together, we arrive at 



_ X 

8=1 



W{e'^^) - iCef + jdwR{w) 
y -iii J 

+ dw R{w) — — - — dw R{w) 



where Fi is the number of massive flavors, F2 is the number of massless flavors, and r is the 
renormalizcd gauge coupling according to 



2mT = 27rzTo + i^Y^^og^ - ^ x {n - \F^ - F^) . (4.27) 



Note the coefficient of the cutoff Q. in this renormalization: The same combination N — 
|Fi — F2 appears in dimensional deconstruction as the coefficient of the 5D Chern-Simons 
term. We are not sure of the physical meaning of this coincidence, but it probably isn't an 
accident. 



36 



Eq. (4.26) gives the effective superpotential in terms of the resolvent R{w), and to re- 
cast it as a function of the gaugino condensates Si,. . . we need to solve eqs. (3.31) for 
the resolvent's parameters. Since the resolvent with n branch cuts has n + 1 independent 
parameters (including the S = '^^Si), determining all of them requires one more equation 
besides n eqs. (3.31). In the matrix model, such n + 1— equation is eq. (3.45), which follows 
from minimizing the free energy with respect to C. The contour integrals in eq. (3.45) are 
similar to a massless flavor's contribution to the To and can be simplified in the same way: 
moving the contours according to fig. (4.22), we arrive at 

^= JdwR{w) - as ^ 0. (4.28) 
-in 

The the effective superpotential W^^{Si, . . . ,5'„) as a function of the gaugino condensates 
follows from combined eqs. (4.28), (3.31), and (4.26). 
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5. Open Questions 



We conclude this paper by discussing open questions raised by our unitary matrix model. 

The first question is specific to the [SU{N)]^ model: what, if anything, is the gauge- 
theory dual of eq. (4.28)? The chiral ring of the quiver (1.3) was studied in detail in [12], and 
none of the anomaly equations there looks remotely like eq. (3.45) or eq. (4.28). Instead, 
the n + 1— equation for the parameters of the R{X) resolvent comes from a completely 
different source, namely the loop equations for the link resolvent T(X) = tr (^ x-nK -Qi ^ ■ 
The analytic properties of the two resolvents R{X) and T{X) require Y{X) {cf. eq. (3.29)) 
to have the same branching points in the X plane as 



HiX - - 4 ((-^)NfA2N.-N,f X n(/if - X) (5.1) 

\ i=i /=i 

where zu^ . . . ,tUj^^ are the Coulomb moduli of the [SU{N)]^ theory; only N^. — I of these 
moduli are independent because of the SU(A'^c) constraint 

J]^. = v^^^. (5.2) 
i=i 

Together, eqs. (5.1-2) impose one constraint on the parameters of the gaugino-condensate 
resolvent R{X)] the remaining parameters are related to the individual condensates Si ac- 
cording to eqs. (3.32), exactly as in the matrix model. 

Unfortunately, the link resolvent T{X) and the Coulomb moduli zu-^ . . . , zuj^^ of the gauge 
theory do not have any matrix duals. Consequently, eqs. (5.1-2) do not make sense on the 
matrix-model side of the duality, just like eq. (4.28) does not make any sense on the gauge- 
theory side. Ideally, these equations could be dual to each other, but we do not have any 
evidence for such duality. In fact, we do not even know whether eqs. (5.1-2) and (4.28) 
are even mutually consistent. This remains an open question we hope to answer as soon as 
possible. 

A bigger open question concerns generalization of our unitary matrix model to other 
quiver theories. We believe that the random unitary matrices can be used to model all kinds 
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of chiral (n, n) bifundamental fields witfi non-zero eigenvalues, but we would like to see how 
this works in different models. Also, it would be interesting to see what exactly goes wrong 
when an eigenvalue happens to vanish in some vacuum state of the gauge theory. We expect 
the gauge-matrix duality to fail for such vacuum, but we are not quite sure, and we certainly 
do not know the specifics (if any) of this failure. 

Finally, we would like to build a matrix model of a chiral (n, ii) bifundamental with 
riy^m. As discussed in the introduction, this calls for complex N x M matrices integrated 
over some variety F C C^^^ which has real dimension NM and satisfies the symmetry 
condition (1.1). More generally, for a theory with several chiral bifundamentals, one may 

use correlated matrices: instead of independent matrices Aii, . . . , Aix each integrated over a 
separate variety F^ {£ = 1, . . . , K), the array {Aii, . . . , Aix) of all the matrices is integrated 
over a combined variety F C C"^ where D = Yle ^i^i. Such combined variety should 
have real dimension D (same as the net dimension of all the F^) and satisfy the generalized 
symmetry condition 

V^f e (net symmetry group), gf : F i-^ F' = F. (5.3) 

We would like to construct such a variety for an interesting quiver theory, and then compare 
the loop equations of the matrix model to the anomaly equations for the field theory's chiral 
ring. 

Acknowledgments: We would like to thank Jacob Sonnenschein and Chethan Krishnan for 
interesting conversations on this subject. V. K. thanks the HEP department of Tel Aviv 
University for hospitality during his multiple visits there. 

This material is based upon work supported by the National Science Foundation under 
Grant Nos. PHY-0071512 and PHY-0455649, and with grant support from the US Navy, 
Office of Naval Research, Grant Nos. N00014-03-1-0639 and N00014-04-1-0336, Quan- 
tum Optics Initiative. Neither funding agency should be held responsible for any errors or 
deficiencies of this paper. 



39 



References 



[1] R. Dijkgraaf and C. Vafa. Matrix models, topological strings and supersymmetric gauge 
theories. Nucl. Phys., B644:3, 2002, hep-th/0206255. 

[2] R. Dijkgraaf and C. Vafa. On geometry and matrix models. Nucl. Phys., B644:21, 2002, 
hep-th/0207106. 

[3] R. Dijkgraaf, M. T. Grisaru, C. S. Lam, and D. Zanon. Perturbative computation of 
gluehall superpotential. Phys. Lett, B573:138, 2003, hep-th/0211017. 

[4] R. Dijkgraaf, S. Gukov, V.A. Kazakov, and C. Vafa. Perturbative analysis of gauged 
matrix models. 2003, hep-th/0210238. 

[5] R. Dijkgraaf and C. Vafa. A perturbative window into non-perturbative physics, hep- 
th/0208048. 

[6] F. Cachazo, M. Douglas, N. Seiberg, and E. Witten. Chiral rings and anomalies in 
supersymmetric gauge theory. JHEP, 12:071, 2002, hep-th/0211170. 

[7] F. Cachazo, N. Seiberg, and E. Witten. Phases ofJ\f = 1 supersymmetric gauge theories. 
JHEP, 02:042, 2003, hep-th/0301006. 

[8] F. Cachazo, N. Seiberg, and E. Witten. Chiral rings and phases of supersymmetric 
gauge theories. JHEP, 04:018, 2003, hep-th/0303207. 

[9] C. I. Lazaroiu. Holomorphic matrix models. JHEP, 05:044, 2003, hep-th/0303008. 

[10] K. Landsteiner, C. I. Lazaroiu, and R. Tatar. Chiral field theories, Konishi anomalies 
and matrix models. JHEP, 02:044, 2004, hep-th/0307182. 

[11] E. Di Napoh and V. Kaplunovsky. Quantum deconstruction of 5D SYM with general 
Uc, Uf and kcs and their moduli space. In preparation. 

[12] E. Di Napoh, V. Kaplunovsky, and J. Sonnenschein. Chiral rings of deconstructive 
[SU{n^)\^ quivers. JHEP, 06:060, 06 2004, hep-th/0406122. 

40 



[13] A. Iqbal and V. Kaplunovsky. Quantum deconstruction of 5D SYM and its moduli space. 
JHEP, 05:013, 2004, hep-th/02 12098. 

[14] N. Seiberg. Adding fundamental matter to chiral rings and anomalies in supersymmetric 
gauge theory. JHEP, 01:061, 2003, hep-th/0212225. 

[15] J. McGreevy. Adding flavors to Dijkgraaf-Vafa. JHEP, 01:047, 03, hep-th/0211009. 



41 



